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ABSTRACT 

The observed size distribution of Kuiper belt objects (KBOs) — small icy and 
rocky solar system bodies orbiting beyond Neptune — is well described by a power 
law at large KBO sizes. However, recent work by Bernstein et al. (2003) indicates 
that the size spectrum breaks and becomes shallower for KBOs smaller than 
about 70 km in size. Here we show that we expect such a break at KBO radius 
~40 km since destructive collisions are frequent for smaller KBOs. Specifically, 
we assume that KBOs are rubble piles with low material strength rather than 
solid monoliths. This gives a power-law slope g ~ 3 where the number N{r) of 
KBOs larger than a size r is given by N{r) oc r^^"^; the break location follows 
from this slope through a self-consistent calculation. The existence of this break, 
the break's location, and the power-law slope we expect below the break are 
consistent with the findings of Bernstein et al. (2003). The agreement with 
observations indicates that KBOs are effectively strengthless rubble piles. 

1. Introduction 

The Kuiper belt, a population of small bodies moving beyond the giant planets, was 
discovered when its first member was found in 1992 (Jewitt & Luu 1993). As of late 2003, 
~800 KBOs have been discovered. Due to KBOs' faintness, however, the size distribution 
of KBOs is well determined observationally only for bodies larger than ~100 km (Trujillo 
et al. 2001, Gladman et al. 1998, Chiang & Brown 1999); their size spectrum is consistent 
with a power law N{r) oc (Bernstein et al. 2003). Numerical studies concluded that 
the differential size spectrum below ~100 km should follow a power law with the slightly 
shallower N oc r~^'^ due to the effects of destructive collisions (Farinella & Davis 1996, Davis 
& Farinella 1997, Kenyon 2002). The results seemed consistent with loose observational 
constraints available on the number of ~20 km and ~2 km KBOs (Cochran et al. 1995, 
Holman & Wisdom 1993). 
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In this context, the deficit in small KBOs observed by Bernstein et al. (2003) was a 
surprise. Using the Advanced Camera for Surveys recently installed on the Hubble Space 
Telescope, they found just 3 KBOs of size ~25-45 km where they expected ~85 such bodies 
based on an extrapolation of the accepted best-fit large-KBO spectrum at the time (Trujillo 
et al. 2001). While this observed decrement of more than an order of magnitude in the 
number of small KBOs clearly indicates a break between 45 and 100 km, the exact break 
position and slope below the break may well be refined by future data on small KBOs. Still, 
the results of Bernstein et al. (2003) are inconsistent with the previously expected small-end 
spectrum N oc or ^ = 3.5, at better than 95% confidence. 

This paper describes a simple self-consistent analytic calculation of the break location 
and the slope below the break. Note that using the N{r) oc size spectrum obtained 
by Bernstein et al. (2003) for large KBOs, we can estimate the size below which coUisions 
between equal size bodies should be frequent to be ~ 1 km — well below the observed break 
location. However, this estimate needs two modifications. First, due to the large velocity 
dispersion in the Kuiper belt, small bodies can shatter much larger objects. Since there are 
more small than large bodies, destructive collisions will occur frequently even for objects 
much larger than 1 km. Second, when collisions are important, they reduce the number of 
small bodies; this in turn decreases the frequency of collisions. Therefore, calculations of 
the effects of coUisions and the size below which collisions are important must be done in a 
self-consistent manner. 



2. Slope of the steady-state distribution 

In order to find the break location self-consistently, we first calculate the power-law 
slope q for a coUisional population of bodies. We assume a group of bodies with isotropic 
velocity dispersion v in which the differential number of bodies of radius r is given by a 
power law dN{r)/dr oc r~^. If we assume that the population is in a steady state and that 
mass is conserved in the collision process, the total mass of bodies destroyed per unit time 
in a logarithmic interval in radius must be independent of size. We can use this condition 
to determine q. 

We assume that the main channel for mass destruction is the shattering of larger 'targets' 
by smaller 'bullets' (Fig. 1). Under this condition, 

pr^ ■ N{r) ■ — ^^p^ ■ ■ V = constant . (1) 

Here p is the internal density of each body and rB{r) is the size of the smallest bullet which, 
on impact, can shatter a target of radius r. V is the volume occupied by all the bodies; 
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their velocity dispersion and therefore their distribution within V are assumed independent 
of body size. When supplemented by a relation between the size of the bullet and that of 
the target, Eq. 1 dictates the size spectrum q. 

This very simple formalism based on conservation of mass captures the essence of 
Dohnanyi's (1969) more elaborate pioneering treatment. Based on laboratory experiments 
which involved solid bodies dominated by material strength, Dohnanyi chose oc r. When 
Tb oc r and N{r) oc r^~'^ are inserted into Eq. 1, we retrieve the q = 7/2 of Dohnanyi and 
several subsequent authors (for example, Williams & Wetherill 1994, Tanaka et al. 1996). 
This slope is much steeper than the best-fit small-end q = 2.3 found by Bernstein et al. 
(2003), who rule out g = 7/2 at better than 95% confidence. 

Indeed, work on the structure of small solar system bodies suggests that many of them 
are gravitationally bound rubble piles rather than solid monoliths. Based on oscillating 
lightcurves of the large KBO (20000) Varuna (radius R > 100 km), Jewitt & Sheppard 
(2002) find that this bodies has density ~ 1 g cm~^ and is therefore unhkely to be sohd. 
However, other effects may also be responsible for the lightcurve shape (see, for example, 
Goldreich et al. 2004). The rotation statistics of much smaller bodies {R ~ 10 km) in 
the more easily observed region between the asteroid belt and the sun also suggest that 
small bodies in the solar system arc rubble piles rather than monoliths. That no small 
asteroids are observed to rotate faster than their breakup speed suggests that those which 
were spun up beyond breakup simply broke apart (Harris 1996); this in tiirn suggests that 
these asteroids have no tensile strength. A study including 26 small near-earth asteroids 
came to similar conclusions about asteroid internal structure (Pravec et al. 1998). The most 
detailed probe available of the structure of small KBOs is research on short-period comets, 
kilometer-sized bodies which are thought to have originated in the Kuiper belt. Work on 
the breakup and impact of comet Shoemaker-Levy 9, thought to be 1-2 km in size, indicates 
that its strength before breakup was ~ 60 dyn cm~^ or less (Asphaug & Benz 1996); a body 
like Shoemaker-Levy 9 would have material strength energy at most about ten times less 
than its gravitational energy. These indications motivate an investigation of the influence of 
negligible material strength on the fragmentation spectrum. 

We might therefore replace the oc r destruction criterion used by Dohnanyi with the 
requirement that the kinetic energy of the bullet be equal to the total gravitational energy 
of the target: 

prlv^ - pr^'vl, (2) 
where v^sc cx: r is the escape velocity from a target of size r, and, again, v is the bodies' 
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Fig. 1. — Schematic of the coUisional cascade: bullets of size rB{R) shatter targets of typical 
size R; these targets break into new targets of size r, which are in turn shattered by bullets 
of size rB{r); and so on. Since mass is conserved in collisions, the mass destruction rate of 
bodies of size R is the mass creation rate of bodies of size r. Steady state then requires that 
the rate of mass destruction be independent of body size. 
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constant velocity dispersion. Then 



rB{r) 




eq 



2/3^5/3 



r, 



eq 




V 



(3) 



Physically, r^q is the size of a body whose escape velocity equals the velocity dispersion of 
the system. When density p ~ 1 g cm~^ and the Kuiper Belt's current velocity dispersion 
f ~ 1 km s~^ arc used, Vcq ~ 10^ km ~ the radius of Pluto. Equivalently, a target of size 
Teg, or roughly Pluto's size, would require a bullet of equal mass to shatter it. Then a body 
smaller than Pluto — that is, virtually any KBO — can be shattered by bullets smaller than 
itself. When we substitute Eq. 3, or essentially the proportionality oc r^/^, into Eq. 1, we 
get the power-law slope 



Recently, O'Brien & Greenberg (2003) extended Dohnanyi's treatment to other destruction 
conditions where scales as an arbitrary power of r; they show that g is a simple function 
of this power so that a range of q values can be obtained from a calculation like Dohnanyi's. 
The simple argument we express in Eq. 1 reproduces their analytic results for q. Eq. 3 is a 
special cases of their general power law which is clearly motivated by energy considerations 
and which leads to the spectrum given by Eq. 4. 



The destruction criterion just discussed neglects any energy loss during the impact 
process. It is then a lower limit on the energy needed to shatter and disperse a given target. 
Indeed, numerical simulations and dimensional analysis of impact events find that in the 
'gravity regime', or target size range where gravity dominates material strength, the impact 
energy needed to shatter a given target lies well above the level indicated by Eq. 3 (Housen & 
Holsapple 1990, Holsapple 1994, Love & Ahrens 1996, Melosh & Ryan 1997, Benz & Asphaug 
1999). Further, the rB{r) scahngs indicated by these studies^ are consistently shallower than 
the one in Eq. 3. With oc r'*, they give 1.37 < a < 1.57 rather than the a — 5/3 in Eq. 3. 

Upon insertion into Eq. 1, the rsir) scalings above give 2.95 < q < 3.11. These values 
indicate a spectrum between the one given by Eq. 4 and Dohnanyi's q — 3.5. This range in 
q is consistent with the best-fit slope q — 2.8 ± 0.6 (95% confidence) derived by Bernstein 
et al. (2003) below the break for the classical Kuiper belt and with the best-fit q — 2.3l5;i 
(bounds of 68% confidence contour) slope they find for the entire Kuiper belt. The value 



Again, we assume a constant velocity dispersion for the coUisional population. 



q = 23/8 



(4) 



3. 



Realistic destruction criteria 
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for the entire belt may be skewed downwards by the scattered Kuiper belt data, which 
include too few faint objects for the scattered belt's small-end slope to be well determined. 
The observed KBO spectrum is thus consistent with the assumption that gravity dominates 
material strength in KBOs of size near rtreak- 

That the simulations give rsir) scalings shallower than that of Eq. 3 implies that the 
energy lost in a catastrophic collision depends on the bullet /target size ratio. As has pre- 
viously been noted (see, for example, Melosh & Ryan 1997), we would expect energy loss 
in the impact of a small bullet on a much larger target. Initially the bullet would transfer 
most of its energy to a volume the size of itself at the impact site; much of this energy would 
escape from the site via a small amount of fast ejecta, though some would propagate through 
the target as a shock. 

Somewhat more quantitatively, we can think of a collision between a very small bullet 
and a large target as a point explosion on the planar surface between a vacuum and a half- 
infinite space filled with matter. The analogous explosion in a uniform infinite material leads 
to the Sedov- Taylor blast wave, a self-similar solution of the first type in which total energy is 
conserved as the spherical shock propagates (Sedov 1946, Taylor 1950). By contrast, a point 
explosion in a half-infinite space is a self-similar solution of the second type (Zel'dovich & 
Raizer 1967); the shock moving into the half-space must lose energy as some of the shocked 
material flows into the vacuum. Also, the nonzero pressure increases the momentum in the 
shock. So as the shock propagates, its velocity should fall off faster than it would have given 
conservation of energy but slower than it would have in the case of momentum conservation. 

We can use these considerations to constrain rsir) scalings for catastrophic collisions. 
We assume that a given target is destroyed if the velocity of the shock wave when it reaches 
the antipode of the impact site exceeds the escape velocity (see, for example, Melosh et al. 
1994). Let the shock velocity decay as f shock where x is the distance traveled by the 

shock. If the energy in the shock were conserved, we would expect /3 = 3/2 from dimensional 
analysis; if the momentum were conserved, we would expect (3 — 3. Then the actual point 
explosion solution must have 3/2 < /3 < 3. The criterion for target destruction is 

/ \ 3-2/3 

fn-W [^) - Gp'r' (5) 

where we have assumed the bullet initially deposits its energy in a volume the size of itself. 
This implies 

1+1//3 7(3 + 1 

rs cx r + , g = . (6) 

The 3/2 < /3 < 3 condition requires 4/3 < a < 5/3 and 23/8 <q < 22/7, both of which are 
satisfied by all of the impact simulation and dimensional analysis results. Holsapple (1994) 
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mentions that energy and momentum conservation should represent hmiting cases for the 
impact process and that laboratory experiments involving impacts into sand, rock, and water 
satisfy those limits. Note that the range in q found in previous studies, 2.95 < q < 3.11, 
spans most of the allowed range for q. This suggests that the catastrophic impact process 
and a depend on more specific details of the collisions such as the equation of state. 

At ~ r there should be no energy loss because the initial energy is deposited in a 
volume of linear size r. Eq. 5 reflects this. Then the Teq expression in Eq. 3 is still valid. 



4. Location of the break 

The above calculation of the size spectrum treats N{r) as constant in time. To maintain 
this steady-state exactly would require the power law to extend to bodies of infinite size, 
which is impossible. To find the range of masses where this assumption holds, we first find 
the size r break of the largest KBO to have experienced a destructive collision after an elapsed 
time r. We equate the timescale for destructive collisions for each KBO of size rbrcak to r 
using Eqs. 1 and 5. To get N{r) we note that bodies of size r > rbrcakj having never collided, 
should be effectively primordial at time r. For their size spectrum we write N{r) = iVor^~^° 
where A^q ~ 4 x 10^^°"^ cm^°~^ from observations (Trujillo et al. 2001). This is equivalent 
to a Kuiper belt with 4 x 10^ bodies larger than 100 km. They are spread over an area 
A ~ 1200 AU^ in the plane of the solar system (Trujillo et al. 2001), so F ~ Av/9, where 
fl = 0.022 yr^^ is the typical orbital angular velocity of the Kuiper belt. With q for the 
slope below the break and, as above, Qq and A^o for the slope and normalization above the 
break, we have 



'"break 



A 



4+9Q-2g 

(7) 



If we set r ^ 4.5 x 10^ yr to be the age of the solar system, take 3/2 < /3 < 3, and use the 
observed qo = 5, we get 

20 km < rbreak < 50 km . (8) 

This is consistent with the observed break position of ~70 km (Bernstein et al. 2003). 
Note that if the system had had the high velocity dispersion assumed above over a time 
considerably shorter than 4.5 Gyr, the break would have occurred at a much smaller KBO 
size. We therefore infer that the Kuiper belt's current excited state has been a long-hved 
phase of at least a few billion years' duration rather than a recent phenomenon. 



The evolution of the total mass and velocity dispersion of the Kuiper belt is a potential 
concern, as the break location depends strongly on both. The mass of the Kuiper belt may 
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have been larger by a factor of ~100 when the solar system was very young (10^ — 10^ years 
old) (see, for example, Kenyon 2002). The collision frequency would have been much higher 
then, so collisions during that period might be expected to have increased the break radius. 
At that time, though, the velocity dispersion of KBO precursors is believed to have been 
just ~1 m/s (see, for example, Goldreich et al. 2002). With this impact velocity, rcq ~ 1 km, 
so only targets of size <1 km can be shattered by bullets smaller than they. As a result, 
coUisional evolution during the early solar system should only have affected bodies of size 
<1 km. The observed break in the spectrum must have been created later. The break 
location could have been affected if there was a sufficiently long period during which both v 
and the Kuiper belt mass were large. 

With ~ 3, the mass contained in bodies of size r <^ r^reak is N{r)pr^ oc r. Since 
the mass destroyed per unit time is independent of body size, the timescale on which coUi- 
sional equilibrium is established is (r/rbreak)'?' <S r. Then the steady-state approximation — 
our assumption that the rate at which N changes is much less than the rate of destruc- 
tive collisions — is self-consistent for r ^ rbrcak- Specifically, as rbreak increases, iV(rbreak) 
decreases — both on a timescale r — and the g ~ 3 spectrum below rbreak follows adiabatically 
(Fig. 2). Our formalism yields the asymptotic spectrum far below rbreak even though the 
system is not in steady state overall, since for r -C rbreak the destruction rate is faster than 
the evolution timescale of the system. Dohnanyi (1969) did not discuss the slow decrease in 
N by which the spectrum differs from a true steady state; he claimed that non-steady-state 
power-law solutions do not exist. Bernstein et al. (2003) conjectured that the disagree- 
ment between their results and Dohnanyi 's calculations might indicate a non-steady-state 
condition in the Kuiper belt. However, the discussion above shows that the fragmentation 
spectrum below rbreak should be unaffected by the system's evolution. 

As for the lower size boundary, the strength hmit derived by Asphaug & Benz (1996) 
implies that material strength dominates gravity at r < 0.3 km. Impact simulations reach 
similar conclusions; they put the threshhold in the 0.1-1 km size range (see, for example. 
Love & Ahrens 1996, Melosh & Ryan 1997, Benz & Asphaug 1999). Below this size thresh- 
hold a different q will apply to an equilibrium collisional population. The changes introduced 
by this effect in the KBO size distribution below ~100 m will affect the size distribution of 
larger bodies through catastrophic collisions. Numerical simulations of collisional popula- 
tions indicate that 'waves' may appear in the size spectrum due to a break in the spectrum 
introduced by a different q (O'Brien & Greenberg 2003). However, the simulations indicate 
that the average slope of the 'wavy' spectrum is not affected; also, the distribution should 
asymptotically approach a ~ 3 spectrum far above this size. 
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Fig. 2. — Temporal evolution of the number of bodies. Here we use g = 3 as a numerical 
example. The sohd hne represents the current KBO size distribution. The dotted line is the 
extrapolation of the large-KBO spectrum to small sizes; we assume this line also represents 
the primordial size spectrum. Dashed lines show the size spectrum at earlier times r = 0.2 
and 1 billion years. Because rbreak increases with time, A'"(rbrcak) decreases with time. The 
evolution of rbreak and A?" (rbreak) is much slower than the rates of coUisional destruction and 
creation below rbreak? so these two rates must be very nearly in balance. Then the steady- 
state approximation is valid in this size range and the spectrum below rbreak follows a g = 3 
power law. 



5. Summciry 

We have derived a self-consistent size spectrum 23/8 < q < 22/7 for a collisional 
population of bodies whose binding energy is dominated by gravity. Wc emphasize that 
this spectrum does not truly represent a steady state; instead, the number density of bodies 
decreases slowly compared to the collision timescale. For the case of the Kuiper Belt, the 
spectrum's small-end power-law slope g ~ 3 and break radius rbreak ^ 40 km agree well with 
those found observationally by Bernstein et al. (2003). Since the power-law slope derived in 
the steady-state approximation depends heavily on the particular criterion for catastrophic 
destruction adopted for the bodies, observations of the KBO size spectrum provide a direct 
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constraint on the bodies' internal structure. The close agreement between this slope and 
break radius and the best-fit values found by Bernstein et al. (2003) suggests that large 
KBOs arc virtually strengthless bodies held together mainly by gravity. Further surveys of 
small KBOs between ~10 and ~70 km in size would better constrain both the exact position 
of the actual break in the size spectrum and the power-law slope below the break. Data of 
this kind would thus confirm or refute our analysis. Such surveys would also allow more 
detailed comparison of the break locations in the classical and scattered KBO populations, 
which should reflect differences in the surface densities and velocity dispersions for those two 
groups. 
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